f I * O I 


|\J 


NASA Technical Memorandum 87294 
AIAA-86-1611 

1 NASA-TM-87294 


Perspectives on Dilution Jet Mixing 


For 


Reference 


J.D. Holdeman 
Lewis Research Center 
Cleveland, Ohio 

and 

R. Srinivasan 

The Garrett Turbine Engine Company 
Phoenix, Arizona 


»ot io 3t iakeii fuom this noou 


LIBRARY mn 


:o 


ijyt) 


EARGLEY research center 
tlgRARY,. NASA 
RAVRTOM, VIRGINIA 


March 1986 





NF01517 




AIAA-86-1 61 1 

Perspectives on Dilution Jet Mixing 

J.D. Holdeman, NASA Lewis Research 
Center, Cleveland, OH; and R. Srinivasan, 
The Garrett Turbine Engine Company, 
Phoenix, AZ 


AIAA/ASME/SAE/ASEE 22nd Joint 
Propulsion Conference 

June 16-18, 1986/Huntsville, Alabama 

For permission to copy or republish, contact the American Institute of Aeronautics and Astronautics 

1633 Broadway, New York, NY 10019 

W'gk- 39436^ 




Perspectives on Dilution Jet Mixing 

by J. D. Holdeman* 

NASA Leu is Research Center , 
Cleveland , Ohio 

and 

R. Srini vasan** 

The Barrett Turbine Engine Company, 
Phoenix, Arizona 


Abstract 

A microcomputer code which displays 3-D 
oblique and 2-D plots of the temperature 
distribution downstream of jets mixing with 
a confined crossflow has been used to 
investigate the effects of varying the 
several independent flow and geometric 
parameters on the mixing. Temperature 
profiles calculated with this empirical 
model are presented to show the effects of 
orifice size and spacing, momentum flux 
ratio, density ratio, variable temperature 
mainstream, flow area convergence, orifice 
aspect ratio, and opposed and axially 
staged rows of jets. 
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= jet-to-mainstream momentum 
flux ratio 
= (DR) (R) 2 

= jet-to-mainstream mass flux ratio 
= (DR) (R) 

= jet-to-mainstream velocity ratio 
= (Vj/U m ) 

= spacing between orifice centers 
= spacing between orifice rows 
= temperature 
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Nomenclature 

= jet-to-mainstream area 
ratio for each row 
= (pi/4) / ( (S/Ho) (Ho/D) 2 ) 


Tj = jet exit temperature 

T m = mainstream temperature 

THETA = temperature difference ratio 

= (T m - T)/(T m - Tj>; Eq. (1) 


= orifice aspect ratio (W/L) 

= (S/Ho) ( ( J) - a > 5 Eq. (3) 

= orifice discharge coefficient 
= orifice diameter 


THETAo 


THETAeb 


= maximum temperature difference 
ratio in a vertical profile; 
see Fig. ,4 

= equilibrium THETA 

= Wj/Wt 


= (D) ( (C„> - n ) 

= jet-to-mainstream density ratio 
= (T m /Tj) 

= duct convergence rate 

= duct height at injection 
location 


THETA m i o'* 


THETAmi „- 


U 


minimum temperature difference 
ratio below the centerline 
(for top injection); see Fig. 4 

minimum temperature difference 
ratio above the centerline 
(for top injection); see Fig. 4 

velocity 


= effective duct height 
= Ho except for opposed rows 
of jets with centerlines 
in-line; see Appendix 


Um = mainstream velocity 

Vj = jet velocity 

w, / w T = jet— to— total mass flow ratio 
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( ( (DR) (J) ) > -"(Co) (Aj / A m ) 


1 + ( ( (DR) (J) ) ■=> (Co) (Aj/A„) 

Wi^ 2 ~ = jet half-width below the 

centerline (for top injection); 
see Fig. 4 

= jet half-width above the 
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centerline (for top injection); 
see Fig. 4 

x = downstream coordinate 

= O at injection plane 

y = cross— stream coordinate 

(radial ) 

= O at injection wall 
= y,; at location of maximum 
THETA in a vertical profile? 
see Fig. 4 

z = lateral coordinate 

(circumferential ) 

= O at centerplane 


Introduction 

The need to design or tailor the 
temperature pattern at the combustor exit 
in gas turbine engines, to maximize engine 
power and life, has motivated several 
investigations on the mixing of jets into a 
ducted crossflow. These include single and 
multiple jets injected into a straight 
duct , ”• o • l ~ v and flow and 

geometric variations typical of most gas 
turbine combustors, namely variable 
temperature mainstream, flow area 
convergence, and double or opposed rows of 
jets, either in-line or 
staggered. 10-10 

The combustor dilution zone 

jet-in-crossflow application is a confined 
mixing problem, with from 10 to 50 percent 
of the total flow entering through the 
dilution jets. The result is that the 
average temperature of the exiting flow may 
differ significantly from that of the 
entering mainstream flow. To control or 
tailor the combustor exit temperature 
pattern it is necessary to be able to 
characterize the exit distribution in terms 
of the upstream flow and geometric 
variables. This requires that the entire 
flow field be either known or modeled. 

From the data of Refs. 1, 11, 14, and 17, 

an empirical model was 

developed 4 " a and 

extended 1 1 • 14 ■ 1 2 for predicting 

the temperature field downstream of single 

and double rows of jets mixing with a 

confined crossflow. It was shown in Refs. 5 

and 18 that this model provides a very good 

predictive capability for the modeled 

variables within the parameter range of the 

generating experiments. 

In addition to the evolution and extension 
of empirical modeling schemes, such as the 
one presented herein, rapid advances have 
been made recently in the capability of 
computational fluid dynamics models and 
their application to complex flows such as 
jet (s> -in-crossf 1 ow. These 

codes offer several advantages over 


empiricism, including the prediction of all 
flowfield quantities (rather than only 
those for which empirical models exist), 
flows outside the range of experiments, or 
flows where empirical assumptions are 
invalid. 

An example of the capability and promise of 
this type of code is given in Ref. 18, 
where temperature field distributions 
calculated using a 3-D elliptic 
viscous-flow code with a standard k-E 
turbulence model 23 are compared 
with measurements from selected cases in 
Refs. 11, 14, and 17, and with 
distributions calculated using the 
empirical model reported therein. The 3— D 
code calculations made in Ref. 18 correctly 
approximated the trends from variation of 
the independent flow and geometric 
variables, but they consistently exhibited 
too little mixing. Although improvements in 
numerics, accuracy, and turbulence models , 
should provide more quantitative 
predictions, there would appear to be a 
continuing need for the empirical model as 
a neai — term design tool. 

An interactive microcomputer program (Apple 
DOS 3.3), based on the model of Ref. 5, was 
used in Ref. 9 to study the effects of 
separately varying the independent flow and 
geometric parameters, and to identify the 
key variables and the relationships among 
them which characterize the mixing. This 
code has been extended to model the effects 
of a variable temperature mainstream, flow 
area convergence, non-circular orifices, 
and double rows of jets, both axially 
staged and opposed. 

Description of the Flowfield 

Figure 1 shows a schematic of the dilution 
jet flowfield for jet injection from the 
top wall. The temperature field results are 
presented in three-dimensional oblique 
views of the temperature difference ratio, 
THETA, where 

(T,„ - T ) 

THETA = (1) 

<T„ - T,> 


A sequence of profile plots of this 
parameter at several locations downstream 
of the jet injection plane, for 
single-sided top injection, is shown in 
Fig. 2- In the 3— D plots the temperature 
distribution is shown in y-z planes normal 
to the main flow direction, x. The 
coordinates y and z are, respectively, 
normal and parallel to the orifice row. 

Note that the jet fluid is identified by 
the larger values of the THETA parameter 
(i.e. THETA = 1 if the local temperature is 
equal to the jet temperature, and THETA = O 
if the local temperature is equal to the 
mainstream temperature). The equilibrium 
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THETA -for any configuration is equal to the 
•fraction of the total flow entering through 
the dilution jets, Wj/w-r. 

The primary independent geometric variables 
are the flow area convergence, dH/dx, 
orifice aspect ratio, AR, the spacing 
between adjacent orifices, S, orifice 
diameter, D, and, for double rows, the 
axial spacing between rows, S*. 

The last three of these are expressed in 
dimensionless form as the the ratio of the 
orifice spacing to duct height, 

S/Ho, the ratio of the duct height 
to orifice diameter, H 0 /D, and the 
ratio of the axial spacing to the duct 
height, S„/H 0 . A 

schematic of the several orifice row 
configurations for which calculations were 
performed are shown in Fig. 3. 


The Empirical Flowfield Model 

The empirical model for the temperature 
field downstream of jets mixing with a 
confined crossflow is based on the 
observation that properly 

non-dimensional ized vertical temperature 
profiles everywhere in the flow field can 
be expressed in the following self— similar 
form 0 


< THETA - THETA™ 


(THETAS - THETA mlr ,~'-) 


-<!n2> (y-y = ) 3 

( 2 ) 

= exp (Wi/ 2 *''-) 2 


where THETA is the local temperature 
difference ratio defined by Eq. <1), and 
THETAc , THETA™ lrl ~, 

THETA™,.™ - , W 1 / 2 *, 

Wi/ 2 - , and y c are 

scaling parameters as shown in Fig. 4. 
Correlations have been developed for each 
of these in terms of the independent 
variables Of, S/D, Ho/D, Z/S, 

X/Ho, S„/Ho, and 

the mainstream temperature, and flow area 
convergence. The correlations used in the 
present version of the empirical model are 
given in Refs. 11, 14, and 17, except that 
the equations describing the effects of 
orifice aspect ratio have been modified. 

The complete set of correlation equations 
is given in the Appendix. 

This empirical model has been implemented 
in a microcomputer program, PROFILES 
version 2.1. The source code was written in 
BASIC, uses the Apple DOS 3.3 operating 
system, and requires > 64K memory. 

This code is a direct extension of that 
used in the study in Ref. 9, and provides, 


as did the previous version, a 
three-dimensional oblique view of the 
temperature field for user supplied values 
of the several independent variables. 

Temperature distributions can be calculated 
at any user — specified downstream location 
Profiles may be be shown as ’cold’ THETA 
distributions, per Eq. (3), where THETA=1 
represents unmixed jet fluid, or as ’hot’ 
THETA distributions where THETA =1 
represents unmixed mainstream fluid. 
Individual profiles from the oblique view 
may also be plotted, either separately or 
overlayed. All profiles shown in this paper 
are ’cold’ THETA distributions? examples of 
’hot’ THETA distributions are shown in Ref. 
9. 

Different perspectives on the mixing can be 
obtained by specifying jet injection from 
the top or bottom duct wall, or with the 
plane between jets, the midplane, rather 
than the plane through the orifice center, 
the centerplane, at the edge in the oblique 
view. Profiles may be shown across a span 
of two or more orifices as desired. 

For all options, the flow and geometric 
variables that must be specified are the 
discharge coefficient, density ratio, 
momentum flux ratio, orifice 
spacing— to-duct height ratio, duct 
height— to— orifice diameter ratio, the axial 
offset between rows, flow area convergence, 
mainstream temperature profile, and orifice 
aspect ratio. Although calculations can be 
performed for most flow and geometric 
conditions of interest, they will be most 
reliable for conditions within the range of 
the experiments on which the correlations 
are based as shown in Table 1. The density 
ratio, momentum flux ratio, orifice 
spacing, and orifice size were the primary 
independent variables. The 
orif ice-to-mainstream area ratio, the 
jet-to-total mass flow split, and the 
parameter coupling the spacing and momentum 
flux ratio, which are derived from the 
primary variables are also given in the 
table. Not all combinations of the primary 
variables in the table were tested; only 
those combinations which are within the 
range given for the derived variables 
represent conditions that are within the 
range of the experiments. 


Results and Discussion 

Figures 5 to 20 show example variations in 
THETA profiles as a function of the several 
independent flow and geometric variables. 
The orifice configurations for which 
profiles were calculated are shown in Fig. 
3, with the profile figure/orifice 
configuration correspondence given in Table 
2. Except for configurations with 
non-circular orifices and profiles spanning 
more than two orifices in a row, the 
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orifice locations are shown schematically 
on the profile figures. Although all of the 
conditions for which profiles are shown in 
this paper do not correspond directly to 
experimental conditions in Refs. 11, 14, 
and 17, many are close, and these are 
identified in Table 2. 

Orifice Size and Spacing. At constant 
total orifice area, changes in orifice size 
and spacing can have a significant 
influence on the THETA profiles. This is 
shown in Fig. 5, where (from left to right) 
jets from closely spaced small orifices 
under-penetrate and remain near the 
injection wall, and jets from widely spaced 
larger orifices over— penetrate and impinge 
on the opposite wall. These profiles are 
all at a downstream distance equal to 
one-half of the height of the duct. 

Momentum Flux Ratio. Figure 6 shows 
a typical increase in jet penetration with 
increasing momentum flux ratio, at the same 
downstream distance. Here the hole diameter 
has been decreased with increasing momentum 
flux ratio to maintain a constant mass flow 
ratio between the jets and mainstream- 

Densitv Ratio. The analyses of the 
experimental data in Ref. 1 suggested that 
the effect of varying the density ratio was 
of second order, for flows with a constant 
momentum flux ratio. This was confirmed, 
over a much broader range of density 
ratios in Ref. 11, and is illustrated by 
the profiles in Fig. 7. In part a), the 
density ratio is varied from .45 (hot jets) 
to 2.2 (hot mainstream) at x/Ho = 

.5 for the same momentum flux ratio and 
orifice geometry as in Fig. 2. The 
similarity of the profiles in Fig. 7 a) is 
in sharp contrast to those in parts b) and 
c) where the density ratio is varied at 
constant values of velocity and mass flux 
ratio respectively. 

Orifice Discharge Coefficient. The 
sixth of the primary independent variables 
which must be specified is the orifice 
discharge coefficient. The profile changes 
which result from variation of this 
parameter from .36 to 1 are shown in Fig. 

8. Note than the THETA values increase in 
response to the greater jet air flow as 
C a increases, but the jet 
penetration and profile shape remain 
si mi lar. 

Coupled Spacing and Momentum Flux 
Ratio. Examination of the experimental 
data in Refs. 1, 11, 14, and 17 revealed 
similar jet penetration over a range of 
momentum flux ratios and orifice spacing 
when these independent variables were 
coupled. For example, low momentum flux 
ratios require large, widely spaced holes, 
whereas smaller closely spaced holes are 
appropriate for high momentum flux 
ratios. 2 > 


As shown in Refs. 13 and 18, jet 
penetration and centerplane profiles are 
similar when the spacing is inversely 
proportional to the square root of the 
momentum flux ratio, i.e. : 

S/Ho = C/ ( (J) -=>) (3) 

For single-side injection, the centerplane 
profiles are approximately centered across 
the duct height and approach an isothermal 
distribution in the minimum downstream 
distance when C=2.5. Values of C in Eq. (3) 
which are a factor of >2 smaller or 
larger than the optimum correspond to 
undei — penetration or over — penetration 
respectively (e.g. Figs. 5 and 6). 

Figure 9 shows profiles for different 
momentum flux ratios, when the orifice 
spacing is adjusted according to Eq (3). 

The profiles in parts a) and b) are at 
downstream distances equal to one— half and 
two duct heights respectively. Clearly, 
similar distributions are obtained when J 
and S/Ho are coupled; however, it 
is also evident that flows with smaller 
momentum flux ratios (larger spacing) need 
greater downstream distances to achieve 
equivalent mixing. 

Orifice Size at Constant Spacing. At 
the desired coupling between orifice 
spacing and momentum flux ratio, the 
profiles remain similar with variation in 
orifice size. The profiles in the top and 
bottom rows in Fig. 10 are at downstream 
distances equal to one-half and two duct 
heights respectively. The orifice diameter 
is doubled going from left to right, 
resulting in a foui — fold increase in the 
ratio of the j et-to-mai nst ream flow. The 
result is that the temperature 
distributions are shifted to higher THETA 
values, consistent with the larger dilution 
air flow, but the shape of the 
distributions remains similar. 

Variable Temperature Mainstream. The 
influence of a non— i sothermal mainstream 
flow on the profiles can be seen by 
examining Fig. 11. In part a), 
under — penetrating jets were chosen to cool 
a hot mainstream near the injection wall; 
in part b) , the jets were positioned to 
cool a center — hot mainstream; and in part 
c>, ovei — penetrating jets were chosen to 
cool a mainstream which is hot adjacent to 
the opposite wall. The empirical profiles 
here were obtained by superimposing the 
upstream profile and the corresponding 
jets— in— an— isothermal mainstream 
distribution. 

This gives a good approximation to the 
experimental data as seen in Ref. 18, but 
it must be realized that this model is onl 
first-order accurate, since with a variabl 
temperature mainstream there can be 
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cross-stream thermal transport due to the 
•flow of mainstream fluid over and around 
the jets (and hence to different y 
locations), and this is not accounted for 
in superimposing the distributions. 

Flow Area Converoence. Figure 12 
shows the effect of flow area convergence 
on the temperature profiles for the same 
orifice geometry and flow condition as in 
Fig. 2. In parts a) and b) the duct 
converges to 1/2 of the injection plane 
height in a downstream distance equal to 
2<H 0 ) and Ho respectively 
(i.e. .25 cm/cm and .5 cm/cm). These 
profiles show increasing jet penetration 
(with respect to the local duct height) and 
slightly increased lateral non-uniformity, 
but these profiles do not support the 
observation 13 that convergence 
improves the mixing. 

Orifice Aspect Ratio. Temperature 
distributions from jets issuing from 
orifice row configurations with bluff and 
streamlined slots of aspect ratio equal to 
2.85 and .35 respectively are shown in Fig. 
13. These configurations are shown in rows 
I and J in Fig. 3, and have the same 
orifice spacing and open area as the 
circulat holes shown in row E. 

The jets from the bluff slots. Fig 13 a), 
penetrate less and are more 
2-dimensional across the orifice 
centerplane compared to the circular holes. 
Fig. 2, whereas the jets from streamlined 
slots. Fig. 13 b) , are highly 
three-dimensional and have deeper 
penetration. At downstream locations both 
the slot configurations and the circular 
holes produce very similar completely mixed 
temperature distributions. 

Slanted Slots. Figure 14 shows 
temperature profiles for the same slot as 
in Fig. 13, oriented at 45 degrees to the 
main flow direction. The orifice area and 
momentum flux ratio for this configuration 
are the same as for the circular holes in 
Fig. 2. 

The experimental data in Ref. 16 for this 
conf iguration shows that the centerplanes 
of the jets shift laterally with increasing 
downstream distance. The empirical model 
profiles in Fig. 14 illustrate this shift 
(c.f. Fig. 2), but the asymmetry of the 
distributions observed in the experimental 
data is not modelled. 

Opposed Rows of In-line Jets. For 
opposed rows of jets of identical orifice 
spacing and diameter, with the orifice 
centerlines in-line, it was shown in Ref. 

13 and 18 that the optimum ratio of orifice 
spacing to duct height is one-half of the 
optimum value for single-side injection at 
the same momentum flux ratio. This is shown 
by the profiles for two orifice 
spacing/momentum flux ratio combinations in 


Fig. 15. The profiles in part b) may be 
compared to the one-side injection profiles 
at the same orifice area and momentum flux 
ratio in Fig. 2. Note that C=1.28 (see Eq. 
(3) ) for these cases, or one-half of the 
optimum value for single-side injection. 

A similar case, but with unequal momentum 
flux ratios on the two sides is shown in 
Fig. 16. Here the total jet flow is the 
same as in Fig. 15 b) , but the plane of 
symmetry, defined by the momentum balance 
betweem the opposing rows, is no longer 
midway across the duct. Note however, that 
this plane of symmetry maintains an optimum 
relationship between the orifice spacing 
and momentum flux ratio on each side. The 
limiting case, i.e. flow from only one 
side, is shown in Fig. 16 b). 

Opposed Rows of Staggered Jets. For 
opposed rows of jets of identical orifice 
spacing and diameter, with the orifice 
centerlines staggered, the optimum ratio of 
orifice spacing to duct height is double 
the optimum value for single-side injection 
at the same momentum flux 
ratio. l3 ' 1B The result of 

implementing this relationship is shown in 
Figs. 17 a) and b). The profiles in part a) 
may be compared to the corresponding 
profiles for one-side injection in part c) 
(this figure is the same as Fig. 2 except 
that the 2 -span of the figure includes four 
orifices rather than two). Note that C=5. 13 
for both parts a) and b) in Fig. 17, 
compared with an optimum value of 2.56 for 
one-side injection. 

It was shown in Ref. 14 and 18 that the 
empirical model does not compare favorably 
with the data in this complex case as the 
fluid dynamic interactions here are not 
amenable to a direct extension of the 
simple Gaussian profile and superposition 
type of modeling appropriate for most of 
the single-side and opposed-jet cases of 
interest. 

Double Rows of Holes Figure 18 shows 
temperature distributions for an orifice 
plate with two in-line rows of jets, for 
three different axial spacings between the 
rows. Profiles are shown at downstream 
distances equal to .25, .5, and 1 duct 
height from the location midway between the 
rows. These profiles may be compared with a 
single row of jets from equally spaced 
holes in Fig. 2. 

It was observed from the experimental 
profiles in Refs. 1 and 16 that the single- 
and double-row configurations have very 
similar temperature distributions, as seen 
also in these profiles. In this case the 
empirical model calculations are derived by 
superimposing the distributions from 
independent calculations of the two rows. 

Figure 19 shows temperature profiles for a 
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double-row con-figuration where the trailing 
row is staggered with respect to the lead 
row. In each row in the figure, profiles 
are given for distances equal to .25, .5 
and 1 duct height downstream from the 
trailing row of holes. As in Fig. IB, 
profiles are shown for varying axial 
offsets, S„/H 0 , from 1 to 

O. Note that the model does not extrapolate 
properly as the offset goes to zero, as is 
evident by comparing Figs. 19 a) and 17 b) 
since the staggered configuration with 
S/H 0 =l becomes a single row with 
S/Ho=.5 when 
S„/Ho=0. 

Temperature profiles are shown in Fig. 20 
for a double-row configuration where one 
row has twice as many orifices as the 
other, but with equal orifice areas in each 
row. In part a), the momentum flux ratio is 
the same for both rows, but is optimum only 
for the lead row. The influence of the 
under — penetration of the trailing row is 
evident even at the farthest downstream 
location. In part b), the momentum flux 
ratio of the lead row is increased 
resulting in over — penetration of these 
jets, but this is compensated for by the 
under-penetration of the trailing row jets 
and a well mixed distribution is approached 
with increasing distance downstream. In 
part c) , the geometric configurations of 
the leading and trailing rows are switched, 
but with the same momentum flux ratios so 
that both rows have the optimum combination 
of spacing and momentum flux ratio. This 
yields a well mixed profile, as expected, 
but the result is about the same with the 
closely spaced row of holes trailing. As 
with the double rows of in-line and 
staggered holes, the empirical calculations 
for this case were obtained by 
superimposing separate calculations for the 
two rows. 

Limitations and Applicability 

Examination of the empirical model results 
in Ref. IB shows that correlation of 
experimental data can provide a good 
predictive capability within the parameter 
range of the generating experiments, 
provided that the experimental results are 
consistent with the assumptions made in the 
empirical model. These models must, 
however, be used with caution, or not at 
all, outside this range. 

Considering the results in Ref. 18 in the 
context of Eq. <3) suggests that in general 
the empirical model provides good 
temperature field predictions for 
single-side injection when 1<C<5. 

Similarly, good predictions are obtained 
for opposed in-line jets provided that 
.5<C<2.5. For opposed rows of staggered 
jets, satisfactory profiles were obtained 
with C approximately equal to 5; but it was 
shown in Ref. 14 that opposed staggered 


jets with lower values of C were best 
modeled using correlations for the opposed 
in-line case. 

The empirical model does not work well for 
impinging flows as the experimental 
temperature distributions are not 
consistent with the assumption of Saussian 
profile similarity in the empirical model. 
The experimental profiles for conditions 
giving optimum mixing in opposed 
staggered-jet configurations are also 
somewhat at variance with the model 
assumptions, and although the agreement 
with the data is satisfactory in these 
cases, this must be considered fortuitious. 

It should also be noted that the form of 
the empirical correlations in the current 
model (and previous versions used in Refs. 
4, 5, 9, 11, 14, and 17) precludes their 
use for semi— confined flows (large 
Ho/D or S/D), single jet flows, or 
flows in which it is known a priori that 
the primary assumptions in the model will 
be invalid. 


Summary of Results 

The temperature profiles herein, calculated 
using at empirical model, for single— and 
double rows of jets injected into a 
confined crossflow, show that: 

1) Variations in momentum flux ratio and 
orifice size and spacing can have a 
significant effect on the resultant 
temperature profiles. 

2) At constant momentum flux ratio 
variations in density ratio have only a 
second-order effect on the profiles. 

3) Similar temperature distributions are 
obtained, independent of orifice diameter, 
when orifice spacing and momentum flux 
ratio are coupled. 

4) The mixing of jets with a variable 
temperature mainstream can be approximated 
by superimposing the 

j ets-i n-an-i sothermal — mai nstream and 
upstream profiles. 

5) Flow area convergence results in 
slightly increased jet penetration and 
cross-stream mixing, but the lateral 
profiles are less uniform that for the 
straight duct case. 

£>) For jets from non-circular orifices that 
are symmetric with respect to the main flow 
direction, the effects of shape are 
significant in the region close to the 
injection plane, but farther downstream 
these geometries yield well mixed 
temperature distributions similar to those 
from equally-spaced, equivalent— area 
circular holes. 
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7) Profiles for jets from 45-degree 
slanted slots shift laterally (z) as a 
function of momentum flux ratio and 

di stance. 

8) For opposed rows of jets, with the 
orifice centerlines in-line, the optimum 
ratio of orifice spacing to duct height is 
one-half of the optimum value for 
single-side injection at the same momentum 
ratio. 

9) For opposed rows of jets, with the 
orifice centerlines staggered, the optimum 
ratio of orifice spacing to duct height is 
double the optimum value for single-side 
injection at the same momentum ratio. 

10) At the same momentum flux ratio, and 
with the same orifice spacing, double rows 
of in-line jets have temperature 
distributions similar to those from a 
single row of equally-spaced, 
equivalent-area circular orifices. 

11) Jets from double rows of orifices of 
different size and spacing, or from double 
rows with orifices staggered, may be 
approximated by superimposing independent 
calculations of the two rows, but caution 
should be exercised using this model for 
very small offsets between the rows. 

The PROFILES v2. 1 code used herein is quite 
capable of calculating flowfields which are 
physically unrealistic, and/or represent 
large extrapolations from the test 
conditions in the data base on which the 
model is based. Some of these conditions 
are flagged in the program with warnings, 
but, of course, the user must always 
exercise judgement is using this or any 
other analytical tool. 
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Table 1. Range of Flow and Geometric Variables 
on which the empirical model is based. 


DR 

J 

S/Ho 

Ho/D 

S„/Ho 

orifice aspect ratio 
variable mainstream THETA 
flow area convergence (dH/dx) 

fij/A m 

Wj/Wt 

C=(S/Ho> (SBR(J) ) 

.5 
, 5 
. 125 
4 

.25 

.36 

O 

O 

.025 

.075 

.5 

to 2.5 
to 105 
to 1 
to 16 
to . 5 
to 2.8 
to .5 
to .5 

to . 1 
to .36 
to lO 

Table 2. 

Correspondence between 
model calculations and 
data in Refs. 11, 14, , 

empirical 
experimental 
and 17. 

Figure 

Orifice 

Test 

Ref. 

(Fig. ) 

2 

E 

I I 1-22 

17 

50 

5-1 

A 

1-5 

11 

20 

5-2 

D 

11-50 

14 

123 

5-3 

G 

1-4 

11 

18 

6-1 

E 

1-1 

11 

12 

6-2 

D 

11-50 

14 

123 

6-3 

B 

1-8 

11 

26 

7a- 1 

E 

1-9 

11 

28 

7a-2 

E 

— 

— 

— 

7a-3 

E 

I I 1-22 

17 

50 

7b— 1 

E 

1-9 

11 

26 

7b— 2 

E 

— 

— 

— 

7b -3 

E 

I I 1—21 

17 

48 

7c— 1 

E 

1-9 

11 

26 

7c -2 

E 

— 

— 

— 

7c -3 

E 

1-1 

11 

12 

8 

E 

— 

— 

— 

9a-l 

G 

1-3 

11 

16 

9a— 2 

D 

11-50 

14 

123 

9a-3 

A 

1-6 

11 

22 

9b— 1 

G 

1-3 

11 

16 

9b— 2 

D 

11-50 

14 

123 

9b— 3 

A 

1-6 

11 

22 

10a— 1 

B 

1-7 

11 

24 

lOa-2 

D 

11-50 

14 

123 

10a— 3 

E 

I I 1-22 

17 

50 

10b-l 

B 

1-7 

11 

24 

10b— 2 

D 

11-50 

14 

123 

10b— 3 

E 

III —22 

17 

50 

1 la-1 

D 

— 





1 la-2 

A 

— 

— 

— 

1 lb-1 

G 

— 

— 

— 

1 lb-2 

D 

— 

— 

— 


B 



Table 2 (con’t) 


Figure 

Orifice 

Test 

Ref. 

(Fig. : 

1 lb-3 

A 




llc-2 

G 

— 

— 

— 

llc-3 

D 

— 

— 

— 

12a 

£ 

— 

— 

— 

12b 

E 

X— 30 

11 

71 

13a 

1 

1 1 1-3 

17 

12 

13b 

J 

XI 1-2 

17 

io 

14 

K 

I I 1—19 

17 

44 

15a 

B/B” 

1 1-7 

14 

26 

15b 

A/A" 

1 1-2 

14 

13 

16a 

A/A- 

11-34 

14 

87 

16b 

A/A- 

1-6 

11 

22 

17a 

H/G- 

11-28 

14 

71 

17b 

C/B- 

11-12 

14 

36 

17c 

F 

— 

— 

— 

18a 

D/D— 

I 11—22 

17 

50 

18b 

D/D— 

1 1 1-6 

17 

18 

18c 

D/D— 

— 

— 

— 

19a 

F 

— 

— 

— 

19b 

H/G— 

1 1 1-9 

17 

24 

19c 

H/G— 

— 

— 

— 

20a 

D/A— 

III-ll 

17 

28 

20b 

D/A— 

III— 15 

17 

36 

20c 

A/D— 

I 11—17 

17 

40 

" top row/bottom 
— front row/rear 

row 

row 
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Appendix — Correlation Equations 


Jet Thermal Centerline Trajectory 
Yc/H-ft, 

= (ai X. 3575) ( J) ■ (S/D) - ** <H. Q /D) — (C d ) - (X/H.*,) - ^ (exp (-b) ) 

where a* = mint ( 1+S/H„d> » 23 

and b= (.091) (X/H..,)=C <H. C ,/S)-(J- “> /3.5) 3 

Centerplane Maximum Temperature Difference Ratio 
THETAc 

= THETA eb +( 1-THETA eb ) C( ai ) (J)-- 3a <C rt ) -“(H.^/D)- 1 (X/H..,) - *• 3^ 
where f = 1 . 15C (S/H.*,) / ( 1+S/H«.d> 3 ■ ® 
and THETA eb = Wj /w T 

Centerplane Minimum Temperature Difference Ratios 

(THETA m i „*) / (THETAc:) 

= 1 - exp(-c*) 

where c~ = (a 3 ) ( - 038) ( J) 1 - *= (S/D) 1 - a (H -c ,/D) -=• a7r (Cc) ■ = s:3 = ! (X/H.c,) 1 - 1 
and a 3 = 1 it (Y = /H„c, + W 1/2 VH.c,) < 1 

= (Hr J /H« q ) 3 - i -' it (Y^/H.^ + W 1/2 */H,. c ,) > 1 

(THETAmir, - ) / (THETAc) 

= 1 — exp ( — c “ ) 

where c~ = (Q) ( 1 . 57) (J) — • 3 (S/D) ~ 1 - •* (H«c,/D) ■ (C d ) • sa (X/H.*,) • v 
and Q = 1 it (Y c /H„ q + W 1/2 */H. q ) < 1 

= expC (.22) (X/H. c ,) = ( (J-=)/5-S/H„ c ,) 3 it (Y c /H. q + W x > 1 

Centerplane Halt-Widths 

= (. 1623) (J) - *«» (S/D) “■ 58,3 (H©/H„«,) -“(Cd) - 1 == 5 (X/H„, q ) • = 

<Wi/a-)/H.«, 

= (.2) ( J) - 13 (S/D) - (H„e,/D)“> 30 (Ho/H wc ,) ■ 55 (Cd) - 033 (X/H_ c , ) ■ ia 
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0-f f-centerpl ane Thermal Trajectory 


Y<= , I /Y c 

= 1 - (4) (Z/S)=(exp(-g>) 

where g = ( . 227) < J) • * 4 *' 7 (S/D) _1 (H.c/D) - (C* ) • = 3 (X/H*,,) - 

0-f f -center pi ane Maximum Temperature Difference Ratio 

(THETAc , x ) /THETAc 
= 1 - (4) (Z/S) = (exp (— d) ) 

where d= ( . 452) ( J) - = 3 <S/D> “ x ■ 03 (H.^/D) • ® 3 (C a > • =» (X/H„c> ■ 0:3 

Of f-centerpl ane Minimum Temperature Difference Ratios 

(THETAmi r. „ **'-) / (THETAi , x ) 

= (THETAmin*'' - ) / (THETAc) 


Of f-centerpl ane Half-widths 
= tWi/2*' ) /H«rc, 


The six scaling parameters, Y c /H. q , THETA C , THETA ml THETA mlr> “, 
and W 1/2 _ /H« q5 are used in Eq. (2) to define the vertical profile at any x,z 
location in the flow. For all except the case of opposed rows of jets with 
centerlines in-line, H.c, in the correlation equations is equal to Ho, 
the height of the duct at the injection location. 


Non— isothermal mainstream 

Double (axially staged) rows of Jets 

Opposed Rows of Jets with Centerlines Staooered 

It was shown in Ref. 18 that these flows can be satisfactorily modeled by 
superimposing independent calculations of the separate elements. This 
is accomplished as follows: 

THETA = C THETA i + THETA* - (2) (THETAi > (THETA*) 3/t 1 - (THETAi) (THETA*) 3 

Note that THETA = THETA ± at any location where THETA* = O (and THETA = THETA 
if THETAi = 0>; and that THETA < 1 (provided that THETAi and THETA* are 
each <_ 1). Also, for the completely mixed case THETAee. is equal to 
the ratio of the total jet flow to the mainstream flow as required. 
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Opposed Rows of Jet with Centerlines In-line 


It was observed in Re-F. 3 that the flow-field downstream of opposed jets 
was similar to that downstream of a single jet injected toward an opposite 
wall at half the distance between the jets. This is confirmed by the 
experimental results in Ref. 14 also. Thus for the symmetric case, 

H„ c , = (Ho) /2. 

In general, these flows can be modeled by calculating an effective duct 
height as proposed in Ref. 12, namely; 

(Hoiqltap = (Ho) ([ (Aj/A m ) (J'°)]top)/ 

( E (Aj /A m ) ( J* =*) 3tap H- C (Aj / A m ) ( J - = ) 3 bat tom) 

and 

fH W q]i;, 0 tt orr , = Ho CH W q]tap 


Flow Area Convergence 

This case is modeled by assuming that the accelerating mainstream will 
act to decrease the effective momentum flux ratio as the flow proceedsX 
downstream, thus: 

J(x) = (J) CH(x)/H 0 3= 

Note that the the trajectory and the jet half-widths are calculated in 
terms of the duct height at the injection location, so must be scaled 
by the inverse of the convergence rate, Ho/H(x), to give profiles 
in terms of the local duct height. 


Orifice Aspect Ratio 

It was observed in Ref. 17 that bluff slots resulted in slightly 
less jet penetration and more two-dimensional profiles than 
circular holes, and that streamlined slots resulted in slightly greater 
jet penetration and more three-dimensional profiles. 

This effect is modeled by using the ratio of the orifice spacing to the 
orifice width, S/W, in lieu of S/D in the correlation equations. 

For rectangular orifices with circular ends; 

S/W = (S/D) Cl + (4/pi ) ( AR -1)3-= if AR > 1 

and 

S/W = (S/D) Cl + (4/pi) (1/AR - l)3-=/AR if AR < 1 

where AR = W/L 

Slanted Slots 

Two effects were noted in the experimental results for slanted slots, 
namely that the centerplanes shifted laterally with increasing 
downstream distance, and the axes of the kidney— shaped temperature contours 
were inclined with respect to the injection direction. The former is modeled 
as a function of momentum flux ratio and downstream distance as: 

dZ/S = sinC (pi/2) (a) 3 

where a = minC 1 , (X/H«„q) ( J/26. 4) ■ ==3 

The rotation effect observed in the experimental data is not modeled. 
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Figure 1. Flow schematic for the 
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with a confined crossflow 
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Figure 2. Temperature profiles with increasing downstream distance 
(S/Ho = .5; Ho/D = 4; J = 26.4; DR = 2.2; Cd = .64). 
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Figure 3. Omif ice Configurations 
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Figure 4. Typical vertical temperature 
profile showing scaling 
parameters in empirical model 
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Figure 5. Temperature profiles with variation in orifice size and spacing at constant 
total area ( J = 26.4; X/Ho = .5; DR = 2.2; Cd « .64) 
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Figure 6. Temperature profiles with variation in momentum flux ratio 
(S/Ho = .5; X/Ho - .5; DR = 2.2; Cd = .64). 
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Figure 7. Temperature profiles with variation in density ratio 
(S/Ho = .5; Ho/D - 4; X/Ho = .&,Cd - .64). 
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Figure 8. Temperature profiles with variation of orifice discharge coefficient 
(S/Ho « .5; Ho/D = 4; J » 26.4; DR = 2.2; X/Ho = .5) 
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b) X/Ho = 2 

Figure 9. Temperature profiles with coupled orifice spacing & momentum flux ratio; 
C = 2.57 (DR =2.2; Cd = .64). 
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Figure 10. Temperature profiles with variation in orifice size at constant spacing 
(S/Ho = .5; J = 26.4; DR = 2.2; Cd = .64). 
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Figure 11. Temperature profiles with non-isothermal mainstream 
(X/Ho = .5; DR = 2.2; Cd = .64). 
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Figure 12. Temperature profiles with flow area convergence 
(S/Ho = .5; Ho/D = 4; J = 26.4; DR = 2.2; Cd = .64). 
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Figure 13. Temperature profiles for bluff and streamlined slots 
(S/Ho - .5; Ho/D = 4; J = 26.4; DR = 2.2; Cd • .64). 
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Figure 14. Temperature profiles for slanted slots 

(S/Ho = .5; Ho/D = 4; J = 26.4; DR = 2.2; Cd - .64). 
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Figure 15. Temperature profiles for opposing rows of jets with centerlines in-line 
(Ho/D = 8; DR = 2.2; Cd = .64). 
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Figure 16. Temperature profiles for opposing rows of jets with centerlines in-line 

and unequal momentum flux ratios (S/Ho = .25; Ho/D = 8; DR = 2.2; Cd = .64). 
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Figure 17. Temperature profiles for opposing rows of jets with centerlines staggered 
(DR = 2.2; Cd = .64) 
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Figure 18. Temperature profiles for double rows of jets with centerlines in-line 

(S/Ho = .5; Ho/D = 5.66; J = 26.4; DR = 2.2; Cd = .64; downstream distances 
from midway between rows). 



X/Ho - .25 


.5 


1.0 



B .5 1 B .5 1 B .5 1 

B=cm-Ti/nn-Tji B*rm-T3/m-TJ3 B=rm-T)/nn-T.n 


a) Sx/Ho = 0 



B .5 1 B .5 1 B .5 1 

B=rm-T3/nfl-TJ3 B=m-T3/[m-T_n B=rm-T3/rm-T.B 


b) Sx/Ho = .5 



B .5 1 B .5 1 B .5 1 

B=rm--n/rm-Tj} fcnn-n/cm-TH B=mi-TJ/[Tn-Ti] 


c) Sx/Ho = 1 

Figure 19. Temperature profiles for double rows of jets with centerlines staggered 
(S/Ho = 1; Ho/D = 4; J = 26.4; DR = 2.2; Cd = .64; downstream distances 
from trailing row). 
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Figure 20. Temperature profiles for double rows of jets with different orifice size 
and spacing in each row; Sx/Ho = .25 (Row 1: S/Ho = .5; Ho/D = 5.66; 

Row 2: S/Ho = .25; Ho/D = 8; Both: DR « 2.2; Cd - .64). 
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